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The propagation of classical wave in disordered media at the Anderson localization transition 
is studied. Our results show that the classical waves may follow a different scaling behavior from 
that for electrons. For electrons, the effect of weak localization due to interference of recurrent 
scattering paths is limited within a spherical volume because of electron-electron or electron-phonon 
scattering, while for classical waves, it is the sample geometry that determine the amount of recurrent 
scattering paths that contribute. It is found that the weak localization effect is weaker in both 
cubic and slab geometry than in spherical geometry. As a result, the averaged static diffusion 
constant D{L) scales like ln(L)/L in cubic or slab geometry and the corresponding transmission 
follows (r(L)) oc InL/L^. This is in contrast to the behavior of D{L) oc 1/L and {T{L)) oc l/L^ 
obtained previously for electrons or spherical samples. For wave dynamics, we solve the Bethe- 
Salpeter equation in a disordered slab with the recurrent scattering incorporated in a self-consistent 
manner. All of the static and dynamic transport quantities studied are found to follow the scaling 
behavior of D(L). We have also considered position-dependent weak localization effects by using a 
plausible form of position-dependent diffusion constant D{z). The same scaling behavior is found, 
i.e., (r(L)) oclni/L^ 

PACS numbers; 42.25.Dd, 42.25.Bs, 72.15.Rn, 72.20.Ee 
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I. INTRODUCTION 

The Anderson localization transition in three dimen- 
sions occurs when the lofFe-Regel criterion Mq < 1 is 
met, where k is the wave vector and io is the bare 
mean free path In the diffusive regime, it is known 
that the averaged transmission coefficient decays with 
the sample size like {T{L)) ~ ^vi/L. In the localized 
regime, it is also known that the geometrical mean of the 
transmission coefficient falls off exponentially with L, i.e. 
{T{L))g ~ exp(— L/^), where ^ is the localization length. 
At the localization transition, it has been predicted that 
{T{L)) ~ 2]. This l/L^ behavior is obtained 

by considering the reduction of the Boltzmann diffusion 
constant, Dq, due to weak localization (WL) effects in 
a spherical volume of size L^. It arises from the contri- 
butions of all recurrent scattering paths returning to the 
origin inside the volume, i.e., those paths of length longer 
than L do not contribute to the reduction of Do- As a re- 
sult, the renormalized diffusion constant, D{L), becomes 
size-dependent and has the form D{L) ~ Dq£q/L at the 
localization transition. Equivalently, the mean free path 
is renormalized to i ~ ^o/L. This renormalization of 
intrinsic transport parameters gives rise to the scaling 
behavior of {T{L)) cx D{L)/L cx l/L^ In one and 

two dimensions, the recurrent scattering paths give di- 
vergent contributions to the reduction of Dq when L is 
large, i.e., ^Z? cx L in ID and 51? cx In L in 2D . Thus 
all states are believed to be localized in one- and two- 
dimensional random media For electrons, inelastic 
scattering due to electron-electron or electron-phonon in- 
teraction provides a nature cutoff length for the recurrent 
scattering paths. Since the inelastic scattering time is in- 



versely proportional to some power in the temperature, 
i.e., Tin oc T~P, WL leads to a temperature-dependent 
conductivity which decreases like —T~p^^ in ID and 
p In T in 2D as long as the dephasing length Ldep oc ^r,;„ 
is smaller than the sample length L. Such temperature- 
dependent of conductivity has been observed in disor- 
dered metal wires and films P, Q. For classical waves, 
the observation of light localization and the scaling be- 
havior of {T{L)) ~ {Iq/LY have been reported ji^l- 
However, these reports have come under close scrutiny 
because of the presence of absorption in the samples. 

Here we would like to point out that the renormaliza- 
tion of Dq discussed above may not be valid for classical 
waves or for electrons when the dephasing length Ldep is 
larger than the sample size L. In this case, the sample ge- 
ometry determines the cutoff of the recurrent scattering 
paths. For classical waves, samples used for the transmis- 
sion measurements are usually not in spherical geometry 
For example, slab geometry are often adapted in 
optical and ultrasonic measurements Since the 
recurrent scattering paths in a slab are different from 
those in a spherical sample, it is natural to ask whether 
the previously obtained scaling behaviors for electrons 
are actually applicable to classical waves? If not, what 
should be correct scaling behaviors for classical waves? 
The purpose of this work is to address these questions. 

In this work, we study the propagation of classical 
waves in finite size disordered samples at the localiza- 
tion transition or mobility edge under the framework of 
the self-consistent theory of localization The con- 

tributions from all recurrent scattering paths within the 
slabs are calculated in the framework of self-consistent 
theory of localization. We show that the averaged static 
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diffusion constant D{L) is proportional to InL/L at the 
mobility edge for both cubic, cylindrical and slab geome- 
tries, in contrast to to the behavior of D{L) (x 1/L ob- 
tained previously for elecrons or spherical samples 0. 
The corresponding static transmission follows the scal- 
ing {T{L)) cx Ini/L^. For dynamics, we have studied the 
time-dependent wave propagation in disordered slabs by 
using both the Bethe-Salpeter equation and the diffusion 
equation with a frequency-dependent diffusion constant. 
In both equations, the effects due to WL are incorpo- 
rated in a way that renormalizes the mean free path . 
It will be shown that the diffusion equation produces the 
same scaling behavior as that of the B-S equations when 

We have also consider the position-dependent WL ef- 
fects by using a plausible form of position-dependent dif- 
fusion constant D{z). We find that the scaling behavior 
(r(L)) cx InL/L^ holds when L/io is large. Since the 
localization effect studied in this work is a general wave 
phenomena, these new scaling laws are not limited to 
classical waves, but may also apply to electrons if Ldep is 
larger than the sample size L. 



II. THEORY 

A. Scaling behavior of renormalized averaged 
diffusion constant at mobility edge 

Weak localization due to interference of recurrent scat- 
tering paths can be signified by the reduction of diffu- 
sion constant in the frequency domain. By summing all 
the maximally-crossed diagrams ,12] in self-consistent di- 
agrammatic theory, the renormalized diffusion constant 
in a bulk can be written as llOj 



scaling behavior that is different from that in spherical 
geometry. The result can be written as 



G{uj; z = z') = sin^ [(?„(z + z,)\ 

n— 1 



-iuj + Dniql+qf,)' 



(3) 



where g„ = m:/ L, = aL/nlQ is the upper momentum 
cutoff in the z-direction, q|| is the momentum in the x-y 
plane, L = L + 2ze is the effective thickness of the slab 
and Ze — 0.7104£o is the extrapolation length Q. We 
let a = 1 in our calculations. A different choice of a 
will only change the mobility edge, kc, not the scaling 
behaviors of wave transport. Eq. (3) indicates that the 
renormalized diffusion constant is z-dependent. In order 
to simplify our calculations, we take the spatial average 
along the z-axis and replace the factor of sin'^[(7„(z -I- zq)] 
by 1/2. This averaging replaces the position dependent 
diffusion constant by its harmonic mean. The situation 
of z-dependent diffusion constant will be considered in 
Sec. lie. 

Here we define a renormalization factor S^iuijk) that 
renormalize the diffusion constant in a finite size slab of 
thickness L according to 



DL{io,k) = 



vIl{l^, k) 



Do 



l+SL{to,k)' 



(4) 



where SL{to,k) = 2TrvG/k'^ by comparing Eq. (4) with 
Eq. (1), and £L{uj,k) = £o/[l + (5^(0;, fc)] is the renor- 
malized mean free path. By using Eqs. (3) and (4), we 
obtain 



1 1 2ttv ~, 

D[uj,k) Do fc^ 



(1) 



where D{uj,k) = v£{uj, k)/3, Dq = vio/5 is the Boltz- 
mann diffusion coefficient and G is the Green's function 
that satisfies the diffusion equation in the frequency do- 



po V + iLu)Giu;; r, r') = -5(r - r'). 



(2) 



The diagonal term of the Green's function, G{w,r,r'), 
represents the return probability of waves that travel dif- 
fusively in the bulk. For an infinite medium, D{u;) has 
already been studied previously in different dimensions 
El in III III 113 ■ To study the scaling behavior in 
a slab, we first solve for G in a cylinder of length L 
and radius to with open ends in cylindrical geome- 
try. A slab can then be obtained by taking the limit of 
R/ L oo. Later we will explain that either the cylindri- 
cal geometry or the cubic geometry would give the same 



k) 



2fc2L 



E 



2q\\dq\\ 



-t^ + Doiql + q^)' 



(5) 



Eq. (5) is then solved self-consistently by replacing Dq 
with DL{uj,k) 0,113. Physically, the self-consistency of 
£L{uj,k) or DL{uj,k) assures the successive renormaliza- 
tion of the recurrent scattering paths inside the samples. 

The scaling properties of disordered slabs can be ob- 
tained by investigating the the scaling behavior of the 
renormalization factor 5l(0, k) in the static limit. When 
i 3> ^Oi the summation in Eq. (5) can be replaced by an 
integral, leading to 



(5^(0, fc) 



2TTk^DL{k) 



dqi 



2q\\dq\\ 



(6) 



where DL{k) = DL{Q,k) and dq± = Ag„. Each of the 
double integrals in Eq. (6) can be split into two parts, 
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and thus Eq. (6) can be expressed as the sum of four 
terms: 



^L(0,fc) 



where 



2'Kk^DL{k) 



(va+Vb +Vc + Vd), (7) 



Dn 



ea\ 2\n{L/£o) 
L J ln2 + 7r/2 



(13) 



Eq. (13) gives D{L) = DL^kc) oc InL/L at fc = kc, which 
can also be obtained from Eq. (4) by using ^i(0,fcc) oc 
L/ln(L). In the following section, we will see that the 
scaling of D{L) dictates the scaling behaviors of many 
measured static and dynamic transport quantities. When 
k > kc, Eq. (13) gives the following L-dependent static 
diffusion constant for slabs: 



r,A ^ (ln2+-)/4 

?7B ~ -(24/i)[ln(L/4) + l] 

77C -(4/i?' + 7r/i?) 

r]D ^ L/R^ + tt/R. 



(8) 



The mobility edge in a bulk can be obtained by taking 
the limits of L ^ oo and i? ^ oo in Eq. (8). At such 
limit, riB — — Vd = and Eq. (7) becomes 



where Doo{k) = Z?oo(0,fc). By substituting Eq. (9) into 
Eq. (4), we obtain the following expression for Doo{k): 



Do 



{io/L)HL/£o) 
(4/Cs)ln(^./4) 



L<^s 
i>6 



(14) 



where is the saturation thickness beyond which DL(k) 
becomes virtually independent of L. can be estimated 
by requiring the L-independent term equal to the L- 
dependent term in Eq. (13), yielding 



4, (is\ 0.72(fc-fc,) 



(15) 



(9) Eq. (15) gives a scaling behavior of 



\k-kc\lTl\k- 

kc\. The above results are different from those obtained 
previously for electrons or spherical samples. For a spher- 
ical sample of radius L, Eq. (2) gives 



npies. t 



Doojk) 
Do 



1 



27r(Ho)' 



(10) 



For the convenience of discussions, here we set £q as the 
units of length and let k to vary. Since the Anderson 
transition occurs when Dao{kc) = 0, Eq. (10) gives the 

mobility edge fc^^o = 
rewrite Eq. (10) as 



3(ln2+;r/2) 
2ir 



1.039. We can now 



D{k) 
Do 



(11) 



For either a 'cubic-like' sample {R = L) or a slab {R 
oo), r]B dominates and Eq. (7) can be approximately 
written as 



SL{0,k) 



Do 



DL{k) 



21n(£/£o) ' 
L J ln2 + 7r/2 



(12) 



By using the relation Do/D^lk) = 1 + 5z,(0, fc) in Eq. 
(12), it is easy to see that 5l(0, k) oc L/\n.L when k — kc- 
By substituting Eq. (12) into Eq. (4), the static diffusion 
constant D^lk) in a finite slab but with L ^ £q can be 
expressed as 



^L(0,fc) 



in^k^DLik) Jl-1 g2 • 
By substituting Eq. (16) in Eq. (4), we obtain 



(16) 



DLjk) 
Do 



1 - 



TT {keo 



1 - 



L 



(17) 



from which we obtain D(L) = Di^{kc) oc £o/ L at k — k^ 
and when k ^ kc, 



DLjk) 
Do 



io/L 



(18) 



where ^ is the correlation length and is proportional to 

\k-kc\-^ m. 

Here we clearly see that the scaling behavior of an in- 
trinsic transport parameter in slabs can be very different 
from that in spherical systems due to different sets of 
recurrent scattering paths. Eqs. (13) and (17) indicate 
that the reduction of D{L)/Do a.t k — kc due to WL in 
slabs is smaller than that in spherical systems by a factor 
of InL. This can be explained by the difference between 
Eq. (6) and Eq. (16) in volume of integration of the 
diffusion pole in momentum space. Fig. la shows the 
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FIG. 1: The schematic diagrams to show the momentum 
spaces of diffusion for the (a) cyhndrical and (b) spherical 
geometry. The shaded areas represent the regions of momen- 
tum space that allow diffusion. For (a), the momentum space 
is a cylindrical volume and the rectangles shown are the sur- 
faces for volume of evolution about the q± axis. For (b), the 
momentum space is a spherical volume and the circles shown 
are the cross-sections of the concentric spheres about the ori- 
gin. 
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FIG. 2: The function /(L) in D{L) ~ /(L)/L at k = kc is 
plotted versus ln(L/^o) for a cubic of length L and a cylindri- 
cal slab of thickness L and radius R = lOL. 



schematic diagram for the momentum space of the cylin- 
drical system. The rectangles shown here are the surfaces 
for volume of evolution about the q_L axis. The shaded 
region represents the momentum space that allows diffu- 
sion, which is schematically equal to the total volume of 
the whole largest square minus the two grey rectangles 
and the small white rectangle. This can also be seen from 
Eqs. (7) and (8). The term of tja represents the largest 
square with a length equal to the upper cutoff of £q^ in 
Eq. (6), which determines the mobility edge, i.e. the 
value of kc- |?7_b| + |r?c| — \vd\ represents the sum of the 



grey and white rectangles, which are to be excluded from 
diffusion. For 'cubic-like' samples or slabs, the value of 
r]B dominates in Eq. (7) and this is the term which makes 
these samples different from the spherical samples. The 
exclusion of the grey regions implies that diffusion is not 
allowed when the momentum in each direction is smaller 
than the lower momentum cutoffs, i.e. g|| < 1/R and 
q± < 1/L. This restriction effectively decreases the total 
return probability and thus reduces the effect of WL sig- 
nificantly. The same argument also applies to the cubic 
geometry, which also involves the separation of variables 
in momentum space into two or more directions. 

However, there is no such restriction in spherical sys- 
tem. The schematic diagram for momentum space of the 
spherical system is shown in Fig. lb. The circles shown 
are the cross-sections of the concentric spheres about the 
origin. The shaded area shows the region of allowed dif- 
fusion modes, which is represented by the volume inte- 
gral in Eq. (16). In spherical geometry, the isotropy of 
the momentum space retains the most recurrent scatter- 
ing paths and thus has the largest WL effect among all 
geometries in three dimensions. At k = kc, the static dif- 
fusion constant depends only on L and it can be written 
as D{L) ~ f{L)/L. Eq. (17) indicates f{L) is constant 
in spherical system while Eq. (13) indicates f(L) oc Ini 
in cylindrical system. Although we do not have a sim- 
ple analytical expression of D{L) for cubic system, we 
also expect f{L) c>c InL in a cube. In Fig. 2, we plot 
the function f{L) versus \n{L/£o) for a cylindrical slab 
of thickness L with radius R = lOL and a cube of length 
L. In this graph we can see that f{L) cx InL and D{L) 
is indeed proportional to InL/L in both cylindrical and 
cubic systems as expected. 



B. The scaling behavior of wave propagation 
through disordered slabs at mobility edge 

To study the scaling behavior in dynamics of wave 
propagation, we consider a pulsed plane wave normally 
incident on the front surface of a disordered slab of thick- 
ness L at z = 0. We assume that there is no gain or ab- 
sorption in the medium and that scattering is isotropic. 
The physical quantity we are interested here is the en- 
semble averaged intensity {I{t, r)), which can be obtained 
from the Fourier transform of the field-field correlation 
function in frequency, i.e., Co(ti^;r) = (r) (/)^_ (r)) 
and 8, 20] 



-iujt)Cn{uj; r) 



(19) 



where $7^ — Q±uj/2, is the central frequency, uj is the 
modulation frequency and 0n(r) is the wave field at po- 
sition r inside the sample with a frequency fl. Cn(w;r) 
can be obtained by solving for the space-frequency cor- 
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relation function in the following B-S equation, 
Co(^;r,r') = {^^^{r)) (<^jV (r')) 

+ J dn dY2 drs dYi{Gn+ (r, ri)) ( G*^- (r', rg)) 
X [/j2(w;ri,r2;r3,r4)Cn(w;r2,r4), (20) 

where ((/)n(r)) is the coherent source inside the sam- 
ple, and (Go(r, ri)) = -^^7(1:377]^ is the ensemble- 
averaged Green's function that represents the coherent 
part of wave propagation from ri to r 1^ . The complex 
wavevector k — k + describes the ballistic propagation 
inside the disordered slab, where fc = — is the wavevector, 
V is the phase velocity, and i is the scattering mean free 
path, which is determined from the imaginary part of the 
self-energy of (G). In the absence of WL, the bare mean 
free path, £0 , is determined from the single-scattering 
diagram via £0 — 1/na, where n is the density of scatter- 
ers and a is the total scattering cross section. The vertex 
function Uq represents the sum of all irreducible vertices. 
Here we approximate Uq as 

;7n(u;;ri,r2;r3,r4) = 

— [1 + Sl{uj, k)] S{ri - rs) ^(ri - rg) <5(ri - r4X.21) 

to 

The first term in the vertex function with a scattering 
strength 47r/£o represents self-avoiding paths and gener- 
ates all the ladder diag rams that give rise to wave diffu- 
sion when L ^ £0 The second term with a vertex 
strength AttSl (o;. fc) /£o represents WL contribution to the 
vertex function . The presence of this term renormal- 
izes the bare mean free path to a frequency-dependent 
mean free path, i.e., ^L(a;, k) = £n/\l + Sr,(LU. k)]. For flux 
conservation to hold, the Ward Identity |lfl llfil| requires 
that the mean free path £ that appears in (G) should also 
be replaced by the same ^^(a;, fc). Since the second term 
represents recurrent scatterings, it is obtained by sum- 
ming all maximally-crossed diagrams due to weak local- 
ization. The renormalization factor (5i(a;,fc) is obtained 
by solving Eq. (2) with appropriate boundary conditions 
and for the slab geometry considered here, we can use 
Eq. (5) by taking the limit of i? ^ 00. 

By using a pulsed plane-wave excitation, the averaged 
intensity is uniform over the transverse cross-section of 
the slab and Eq. (20) can be expressed as 



6n(w, z) = exp 



exp 



V £l(0) 



1 



47r£L(w,fc) Jo 
where ^l(O) = ^L(0,fc) and 



dz'H{uj,z- z')Cn{oj,z') ,(22) 



p2 + (z — z'Y 



(23) 



Eq. (22) is then numerically solved for Cn{Lu,z). The 
transmitted intensity (/(t, L)) is calculated from the 
Fourier transform of Cn{i^,L) as in Eq. (19) and 
Gn(0, L) gives the static transmitted intensity {I{L)). 

In the second method to obtain {I{t,r)), we solve the 
diffusion equation in the frequency domain with a fre- 
quency dependent diffusion constant, DL{uj,k). The WL 
effects are incorporated through DL{uj,k) according to 
Eqs. (4) and (5) with R 00. The solution to the 
diffusion equation in a slab takes the form 



2 ^ sin[g„(z, 

^ n=l 



Zp)] sin[q„(z -|- Ze)] 



iuj + Dl{uj, k)ql 



(24) 



where q„ = nil / L is the transverse momentum, L = 
L + 2ze is the effective length, Ze — 0.71£l(0, fc) is the 
extrapolation length and Zp ~ ^l(0, fc) is the penetration 
length. Here we use the renormalized mean free path 
£l{Q, fc) — £q/{1 + 5l(0, fc)) in the evaluation of the ex- 
trapolation length Ze and the penetration length 'Zp. This 
replacement is consistent with the replacement of £0 by 
£L{i^,k) in both the averaged green function (G) and ver- 
tex function Uq in the B-S equation as required by the 
Ward identity. It should be mentioned that the dynamic 
diffusion constant D{uj) has been studied for electrons 
near the mobility edge In these studies, a be- 

havior of D(lj) cx u}^^^ was found at mobility edge for an 
unbounded medium. 

The static transmitted intensity {I{L)) can be ob- 
tained from Eq. (24) by setting a; = 0, yielding 



im) = 

The transmission {T{L)) 
the form 



Ze ~t~ 



L + 2zeDL{ky ^ ' 

= -DL{k)±{I{z))l^L and has 



{TiL)) = 



2z. 



(26) 



At fc = fcc, Eqs. (25) and (26) give the same scaling 
behavior of InL/L^ for both (/(L)) and (T(L)). It is 
worth to note that it is the replacements of Zp ~ by 
Zp ~ £i(0, fc) and Ze 0.71^ by Ze ~ 0.71^^(0, fc) that 
change the scaling of {I{L)) and {T{L)) from to 
lnL/L2. 



C. Discussion on the position-dependent diffusion 
constant 



H{uj,z-z') = 



dp^ 



In an open system, the WL effects should vary in space 
as the probability of returning to each point inside the 



6 



0.8 



0.6 



Q 0.4 



0.2 



relation 1211 



finite slab (L=10001J 

semi— infinite medium 



50 100 150 200 250 



FIG. 3: D{z)/D(Q) for a finite slab witii L = fOOO^ (solid 
line) and the semi-infinite medium (dashed line), which are 
obtained from Eqs. (32) and (31) respectively, are plotted as 
functions of z/£o. 



sample can be different. Thus in Eq. (1) the renormal- 
ized diffusion constant D can also be position dependent. 
In arriving Eq. (12) for slab geometry, we have simplified 
the calculation by taking the spatial average of G along 
the z-axis and assumed that the diffusion constant is in- 
dependent of z. Effectively, this simplification replaces 
D{z) by its harmonic mean, i.e., D{L) = {1/D{z))^^ as 
can be seen from Eq. (1). A complete theory requires 
self-consistent solutions of both D{z) and G(r,r'), from 
which one can obtain {T{L)). In the case of slab geome- 
try and static limit, Eq. (2) is replace d by the following 
position-dependent diffusion equation |2l| . i.e.. 



d_ 

dz 



D{z) 



(iG(q|| ; z, z') 
dz 



where 



G(r,r') = 



(27r)2 



^Diz)qlGi<if,z,z') = -Siz-z'), 

(27) 

dqy exp(iq|| • p)G'(q|| ; z, z'). (28) 



In an operator form, we can write Eq. (27) as 
£{G(q||; z, z')} = —S{z — z'), where the self-adjoint hnear 
operator £ is given by 



{T{L)) - (z, + zp) 



22. 



dz ^ 
Di^) 



(30) 



where the extrapolation length z. and the penetration 
length Zp are determined by the diffusion constant at 
sample boundary |2ll |. By repeating the same calcula- 
tion at different L's, one can obtain the scaling behav- 
ior of (T(L)). For L > io, Eq. (30) can be written 



as {T{L)) ~ (z, + Zp) 



Do 
D{z) 



z^+Zp D{L) 



By 



- Do L 

comparing Eq. (30) with Eq. (26), we can see that 
the scaling of {T{L)) for both position-dependent and 
position-independent diffusion constant are determined 
by the scaling of D{L). However, there is a subtle differ- 
ence between the two situations. In Eq. (26) the scaling 
of {T{L)) is affected by the scaling of Zg and Zp through 
the scaling of II or D{L). However, in Eq. (30) the scal- 
ing of (r(L)) arises directly from D{L)/L. Since {T{L)) 
is dominated by the small values of D{z) deep inside the 
sample, an accurate numerical calculation for {T{L)) is 
difficult when L is large. In this work, we do not intend 
to solve this problem self-consistently. Instead we would 
like to propose a plausible form oiD{z) and show that the 
corresponding transmission is close to the self-consistent 
solution and behaves like {T{L)) cx h\L/lJ^ when L is 
large. 

For a semi-infinite medium, van Tiggelen et. al. 
have performed the self-consistent calculations discussed 
above and suggest an analytical form of D{z) at mobility 
edge lij, i.e.. 



Doo(z) 



(31) 



where I?oo(0) is the diffusion constant at the boundary 
of a semi-infinite medium. Eq. (31) shows Doo{z) de- 
creases like 1/z from its value at the boundary when z 
is moving into the semi-infinite medium. In the absence 
of internal reflection, they find D^{0)/ Dq — 0.642 and 
Cc/^o = 1-5. For a finite slab of thickness L, based on Eq. 
(31), they have also suggested that Dl{z) — Doo{zl), 
where zl = -j — — z| jiJi, from which they found 
{T{L)) oc 1/L^. In the case of finite slabs, we expect 
that Dl{z) should decrease slower than 1/z as z moves 
well inside the sample due to reduced WL effects in the 
presence of other boundary. By taking this into account, 
we propose here the following modified form for Dl{z): 



dz 



D{z)ql 



(29) 



For a rigorous approach, one should solve Eq. (1) and 
Eq. (27) simultaneously with appropriate boundary con- 
ditions at two surfaces for each fixed sample thickness L. 
The transmission coefficient is then obtained from the 



Dl{z) = Z?oo(0) 



{zL/W\n{zL/io + w)-l) 
(ln(zi/4 + w)y 



(32) 

with L ^ w ^ S^c- The above Dl{z) is symmetric with 
respect to the central plane at L/2. Its value decreases 
monotonically from the boundary to its minimum at the 
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FIG. 4: (T<^^(L)> obtained from the first iteration of Eq. 
(33) (solid line) and {T''°\L)) obtained directly from Dl{z) 
in Eq. (32) (dashed line) are plotted against L/lo- In the 
inset, {r(^^)"^-(lnL) (solid line) and (T^"')-^ ■ (In L) (dashed 
line) are plotted versus L/£o in log-log scale. The slope of 
their linear fitting lines are 1.95 and 1.98 respectively. Also 
plotted in the inset is (t''^')"^ ■ (In L) obtained from Eq. (26) 
(dotted line) with a linear fitting line of slope 1.97. 
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FIG. 5: The quantity Eo{ci\\)/ D{L) for L/fo = 500 and 1000 
are plotted against gy^o- In the inset, the curves near £q ~ 
for L/£o = 500 (solid line) and 1000 (dashed line) are ampli- 
fied by a factor of 10'' to show the y-intercepts. The fit to the 
curves suggest that £o(q||) - Eo{0) + D{L)Pqf with fJ ~ 0.55. 
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center. In this work, we choose w ~ 10.65 in Eq. (32) to 
match the behavior of Eq. (31) near the sample bound- 
aries with fc/^o — 1-5. In order to make comparison 
between Eqs. (31) and (32), In Fig. 3 we plot both 
Dl{z)/ Doo{0) of Eq. (32) for the case of a finite slab 
with L/io = 1000 (solid line) and Dodz)/ Doo{0) of Eq. 
(31) for semi-infinite medium (dashed line). From Fig. 3, 
it is easy to see Dl{z) decays like Doo{z) near the bound- 
ary, but decreases in a slower rate than 1/ z away from the 
boundary. Since the value of {T{L)) is dominated by the 
small values of Dl{z) inside the sample as it can be seen 
from Eq. (30), a scaling behavior which is different from 
{T{L)) (X l/i^ is expected. By substituting Eq. (32) 
into Eq. (30) and setting Zp = 3DooiO)/v = 0.642^ and 
Ze = 0.71zp ji^l, we obtain the transmission, (r°(i)), 
which is shown by the dashed curve in Fig. 4. In order 
to show more clearly its scaling behavior, we replot the 
function {{T'^ {L)))~^ InL in the inset of Fig. 4 using log- 
log scale. A linear line of slope 1.98 clearly shows the 
scaling relation {T{L)) cx InL/L^. 

In order to test the self-consistency of this result, we 
substitute Eq. (32) into Eq. (29) and solve for the eigen- 
value problem of £{J7„(q|| , z)} = i^nfc[ji ){C^„(q|| , z)} by 
using the method described in Ref. |22|. The diagonal 
Green's function G(r, r) of Eq. (27) can be obtained by 
the eigenf unction expansion: 



§^10-^ 
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FIG. 6: The values of Eo{0)/D{L) are plotted against L/£o in 
log-log scale. The slope and the y-intercept of the linear fitting 
are -1.99 and 2.59 respectively, thus£'o(0) ~ 1.35D{L){Tv/L f . 



(2^)^ 



ji=0 



\ Uni<l\\,z) 

£^n(qii) 



(33) 



In Sec. IIA, we have shown that the InL factor found 
in rjB of Eq. (8) arises from the exclusion of phase space 
volume bounded by the lowest order mode q± — tt/L 
when qy ~ as shown in the grey area of Fig. 1. Here 
we would like to show that this volume exclusion effect 
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holds even when a position-dependent D{z) is consid- 
ered. Since we are interested in the scahng of {T{L)), 
only the spatial average of G{z, z) is relevant as can be 
seen from Eq. (30). The normalization of eigenfunctions 
requires |J7„(q||)p = 1/L for all n and q||. We can see 
that the forbidden volume in phase space is determined 
by £'„(q||). For a fixed q||, the lowest order eigenval- 
ues for different slab thickness L divided by D{L), i.e. 
Eq{ci\\)/D{L), are plotted versus qjiQ in Fig. 5. The y- 
intercepts for each curves are also shown in the inset of 
Fig. 5. The linear fits to the curves for each L suggest 
that Eo{q\\) ~ Eo{0) + D{L)Pq^^. For sufficiently large L, 
P ~ 0.55. By using the y-intercepts of the fitted fines in 
Fig. 5, Eo{0)/D{L) are also plotted against L/io in log- 
log scale in Fig 6, which are well fitted by the formula 
^0(0) ~ 1.35l)(i)(7r/L)2. The combined results indi- 
cates that Eo{q\i) ~ D{L)[1.35{tt/ L)'^ + /3g^] and the ex- 
clusion of phase space volume is retained for finite L. By 
using this result in Eq. (33), we would expect to obtain 
the same scaling behavior for the transmission. In order 
to confirm our assertion, we first calculate the averaged 
diffusion constant, D'^\l) from Eq. (1) by using the 

results shown in Figs. 5 and 6. By substituting D^^\l) 
into Eq. (30) we obtain (T(i)(L)). This result is plotted 
as the solid curve in Fig. 4. The excellent agreement 
between {T^°\L)) and (r(i)(L)) when L > lOO^o indi- 
cates that Eq. (32) is close to the self-consistent solution 
when L/£o is large. This result strongly indicates that the 
scaling behavior of D{L) ex In L/L or {T{L)) cx Ini/L^ 
found in the previous section holds even when a position- 
dependent diffusion constant is considered. In fact, we 
will show in Sec. Ill that the transmission coefficient 
shown in Fig. 4 agrees very well with that obtained from 
an averaged diffusion constant D{L) given in Eq. (14) 
through the use of Eq. (26) , which is shown by the dotted 
line in the inset of Fig. 4. In the next section, we present 
the numerical results of wave propagation through dis- 
ordered slabs at mobility edge based on the averaged 
frequency-dependent diffusion constant Z?i(w, fc) shown 
in Eqs. (4) and (5). 



III. NUMERICAL RESULTS AND DISCUSSIONS 

Before presenting the dynamic results, we would like to 
show that Eq. (25) of the diffusion approximation (DA) 
is capable of producing the results of the B-S equation 
when L ^ io. We first consider a case in the diffusive 
regime with k > fee- In Fig. 7, we plot the recipro- 
cal of static transmitted intensity (/(i))~^ at k£o = 8 
against the dimensionless slab thickness L/£o in log- log 
scale. Both the results from the B-S equation and the 
DA are shown. It can be seen that the results from the 
DA agrees well with those from the B-S equation apart 
from a small constant shift. The dashed line with a slope 




FIG. 7: The reciprocal of static transmitted intensity 
(/(L))~^ for kio = 8 is plotted as a function oi L/£o in log-log 
scale. Both B-S results and diffusion approximation (DA) re- 
sults are shown. The dashed line with slope of 1 is plotted to 
shown {I{L)} oc 1/L. The scaled reciprocal static transmitted 
intensity {I{L))~^ ■ (InL) for fc = fcc is also plotted with the 
logarithmic scale on the right. The solid line with slope of 2 
shows that {I{L)) oc In L/L^ a,t k — kc- 



of 1 is plotted to show that {I{L)) cx 1/L for both the B- 
S and the DA results. At the mobility edge, i.e. k = kc, 
we would expect {I{L)) cx \nL/L^. On the right scale of 
Fig. 7, we plot the calculated result of {I{L))^^ ■ (InL). 
The solid line with a slope of 2 is also plotted to show 
that (/(L)) calculated from both the B-S equation and 
the DA indeed give the scaling behavior of In L/L^. We 
have also calculated the transmission coefficient using Eq. 
(26). These results are denoted as (T'^"^)) and plotted as 
the dotted line in the inset of Fig. 4. The excellent agree- 
ment between (T^^') and (Tt"') or (T(i)) supports the 
use of an averaged diffusion constant in the transmission 
calculations. 

In Fig. 8, we plot {I{t)) = {I{t,L)) for L/£o = 20 at 
k = kc obtained from the B-S equation (in solid line) and 
the DA (in dashed line) as functions of time t/rn, where 
Tfj = /tt^D{L) is the renormalized diffusion time. Our 
results show that they agree with each other for a rather 
large range of t/ru, indicating the validity of Eq. (24) 
when £o. In the same graph, we also plot the result 
obtained from Eq. (12) of Ref. [l^, in which a time- 
dependent diffusion constant D{t) = Do{£o/tvy^^ is used 
in the time-dependent diffusion equation. Since such a 
local scaling approach does not consider the retardation 
effect of the recurrent scattering paths, it overestimates 
the WL effects and, therefore, produces a smaller {I{t)) 
and with a slower decay rate as shown in Fig. 8. Such 
approach has also been used in the study of coherent 
backscattering [2^ and absorbing media [23| near mobil- 
ity edge. In our theory, it is the frequency-dependence 
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FIG. 8: The logarithmic time-resolved transmitted intensities 
ln(J(t)) at fc = fcc for L/io = 20 obtained from the B-S equa- 
tion (in solid line) and diffusion approximation (in dashed 
line) are plotted as functions of t/rR. The numerical result 
from the local diffusion approximation (in dot-dashed line) 
with the same parameters is also plotted for comparison. 



of the factor 6l in Eq. (5) that makes the reduction 
of intrinsic diffusion constant causal in time. It is also 
interesting to point out that, unlike a pure diffusion pro- 
cess, the decay of {I{t)) shown in Fig. 8 is not a sim- 
ple exponential decay. The slowdown of decay rate in 
time is a result of increasing WL effect contributed by 
the presence of longer recurrent scattering paths. Such 
non-exponential decay has also been reported first for 
electronic systems |26|, |27| and recently observed in the 
microwave experiments in a nominally diffusive region 
|28| . Here the cause of the non-exponential decay is also 
due to WL effects, but in a quasi-lD geometry [llLl27ll2g | 
The non-exponential decay shown in Fig. 8 becomes 
more transparent in the time-dependent diffusion coeffi- 
cient D{t), which is defined via the relation D{t)/Do = 
—td d\n{I{t, L))/dt, for t > tr, where td = th/ (1 -I- 6l) 
is the diffusion time [2^. In Fig. 9 we plot D(t)/Do at 
k = kc, obtained from the DA, for various sample thick- 
nesses L/£o as functions of t/rji. It is easy to see that 
D{t) I Do is not a constant in time and its maximum value 
decreases as L increases. This is in contrast to the case 
of the diffusive regime, i.e. when k > kc, in which the 
WL effect is weak and the change of slope in D{t) is very 
small |3^. 

Since the static diffusion constant D{L) has a scaling 
that is different from that of the spherical samples, it 
is also interesting to investigate the scaling of the time- 
dependent diffusion coefficient D{t). Here we are inter- 
ested in D{t) in the long time limit, i.e. D{t 00), 
because the WL contribution of the long recurrent scat- 
tering paths is expected to saturate eventually. In the 



FIG. 9: The time-dependent diffusion coefficient D{t)/Do at 
k = kc for L/lo = 20, 30, 40, 50, 80, 100 and 500, obtained 
from the diffusion approximation, are plotted as functions of 
time t/TR. 



inset of Fig. 10. we show the fitting of the long time 
tail of D{t)/Do (in solid line) against the function a + b/t 
in dashed line for the case of L — 20£o aX k — kc- We 
perform the similar fitting to each of the curves in Fig. 9. 
From these fittings, we obtain D(t ^ 00) = a for differ- 
ent slab thickness L/£q and they are used to plot D{t ^ 
oo)/DolnL as a function of L/io in log-log scale in Fig. 
10. The solid line is the fit to the curve with a slope of 
—0.99, which suggests that D{t 00) /Do oc \nL/L and 
is consistent with the scaling of D{L) given by Eq. (13). 

Figs. 9 and 10 show that the InL factor appears when 
waves see the boundary of the sample, i.e., t ^ tr. When 
t < tr, waves have not reached the output surface and 
the InL factor should not appear. In fact, we find D{uj) 
behaves like w^^'^ when ui > I/tr. This is consistent with 
the previous studies 0,^3- When u < 1/tr, we find 
D{lu) — D{L) (X ijj^ , which represents the saturation of 
WL effects, or equivalently, the saturation of D(i) when 
i/^i? ^ 1 as shown in the inset of Fig. 10. 

Besides the static intensity and the diffusion coeffi- 
cient, the time of peak intensity Tpcak, i-e. the time 
it takes waves to diffuse across the sample, also follows 
the scaling of <5l(0, /c). Tpcak can be obtained from the 
peak position of the time-resolved transmitted intensity 
(I(t)). In standard diffusion theory, Tpcak — L'^/iDo 
[19j. If weak localization in included, the Boltzmann dif- 
fusion coefficient Do should be replaced by D^ik) and 

thus Tpcak ^ [1+ 5l)L^/?,Do = (1 + 5L)L''/vio- Tpcak 

versus L/io at kio = 8 in log-log scale is plotted in Fig. 
11 with the logarithmic scale shown on the right. Both 
B-S and DA results are shown and they agrees with each 
other except for small sample thickness. The slope of 
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FIG. 10: The scaled long time diffusion coefficients D{t 
oo)/ Doln L at k = kc, obtained from the fitted values of a 
from the curves in Fig. 5, arc plotted as a function of L/£o 
in log-log scale. The slope of the fitted line is -0.99. In the 
inset, the dashed line shows the fitting of D{t — > oo)/Do (in 
solid line) by the function a + b/t for the case of L = 20£o- 
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FIG. 12: The position of peak intensity Zpeak{L) / L obtained 
from the B-S equation is plotted against \ftjTR. The dashed 
line is the fitting to the linear region, which has a slope of 
0.33. 
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FIG. 11: The time of peak intensity Tpeak{L) at kio = 8 
is plotted against L/£o in log-log scale, with the logarithmic 
scale shown on the right. Both B-S results and diffusion ap- 
proximation (DA) results are shown. The slope of the fitted 
line (dashed line) is 1.99. The scaled time of peak intensity, 
Tpeak InL, at k = kc versus L/£o is also plotted with the log- 
arithmic scale shown on the left. The slope of the fitted line 
(solid line) is 3.01. 

the dashed line is 1.99, which confirms that Tpeafe oc 
at weak scattering Umit. The scaled peak time of inten- 
sity, Tpeak^i^L, at fc = kc versus L/£o is also plotted on 
the same graph with the logarithmic scale shown on the 
left. The slope of the fitted line is 3.01, which indicates 
that Tpeak oc L^/\nL and is consistent with the scaling 
of (5l(0,A:c) cx L/lnL. 



We have also studied the position of the peak inten- 
sity Zpeak obtained from B-S equation, which represents 
the collective transport of diffusive waves in the sample. 
If weak localization is present, the recurrent scattering 
paths enhance the backscattering of waves and thus de- 
lay the propagation of the peak intensity. Zpeak{L) / L at 
fc = fcc for L/£o = 20 and 50 are plotted as functions 
of sJiIt-r in Fig. 12. The fitted line with a slope of 
0.33 is added to highlight the region of diffusive trans- 
port. This suggests that Zpeak/L ~ 0.33\/t/TH and thus 
Zpeak ^ \/D{L)t for t < tr, which is consistent with 
the diffusion theory. When t > tr, the curves start to 
deviate from the dashed line because the diffusive waves 
have reached the open end. As we can see in Fig. 9, the 
value of Zpeak/L is always less than 0.5 and this implies 
the wave interference under strong scattering can stop 
the peak intensity from approaching the central plane of 
the slab. 



IV. CONCLUSIONS 

In summary, the scaling behavior of wave transport at 
the Anderson transition, i.e. k = kc has been studied. 
We found that both the static and dynamic transport 
properties at k = kc follow the scaling of the averaged 
static diffusion constant D{L). D{L) is found to scale 
as InL/i in cubic, cylindrical samples or slabs, in con- 
trast to the scaling of D[L) oc 1/L found previously for 
electrons or spherical samples. The corresponding static 
transmission {T{L)) scales like InL/L^, in contrast to 
the l/Li^ behavior found previously. Our results indi- 
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cates that the weak localization effects in other geometry 
are in general weaker than that of the spherical system 
by a factor of 1/lni. This factor arises from the exis- 
tence of larger volume exclusion in the phase space of 
allowed diffusion modes in non-spherical samples. For 
dynamic transport, we solved both the Bethe-Salpeter 
equation and the diffusion equation with weak localiza- 
tion included self-consistently. The numerical results cal- 
culated by the two methods agrees qualitatively with 
each other and they are found to produce the same scal- 
ing behavior when L ^ £o- The scaling for the long 
time diffusion coefficient and the time of peak inten- 
sity are D{t oo)/Dq cx InL/L and Tp^ak oc L^/lni 
respectively, which are consistent with the scaling of 
D{L) in slabs. In addition, the position of peak inten- 
sity in a slab Zpeak/L is found to scale as ^Jt/rn when 
t < Tji, which is also consistent with the diffusion theory 
of D{L). We have also studied the position-dependent 
weak localization effects by using a plausible form of 
position-dependent diffusion constant D{z). The same 
scaling behavior is obtained for the transmission, i.e., 
(r(L)) cx lnL/L2. 
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Hong Kong RGC Grant No. HKUST 6058/02P. 
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